Active control of qubit-qubit entanglement evolution 
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In this work, we propose a scheme to design the time evolution of the entropy of entanglement between two 
qubits. It is shown an explicit accurate solution for the inverse problem of determining the time dependence 
of the coupling constant from a user-defined dynamical entanglement function. Such an active control of en- 
tanglement can be implemented in many different physical implementations of coupled qubits, and we briefly 
comment on the use of interacting flux qubits. 
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I. INTRODUCTION 

In quantum information science, entanglement is com- 
monly recognized as a resource that may be used to perform 
novel information processing tasks 1 1]. A famous example of 
such a task is quantum teleportation 12|] . Two parties sharing a 
pure maximally entangled pair of qubits can perform the per- 
fect transmission of the state of a third qubit using just local 
operations and classical communication. This would be im- 
possible with the use of classical communication alone. This 
view of entanglement as a resource has triggered the investi- 
gation of measures of entanglement. Complete and didactic 
reviews of the main entanglement measures and their funda- 
mental properties can be found in Jl], . 

The common approach to generate entanglement between 
qubits or more complicated quantum systems consists of let- 
ting the system evolve under some time-independent interac- 
tion Hamiltonian. In general, the total time of evolution and 
the initial preparation of the state of the system are the in- 
gredients varied to achieve the desired final entangled state. 
In this passive approach, the system entanglement follows a 
unique time evolution once the initial state is fixed. On the 
other hand, the use of time-dependent Hamiltonians could al- 
low for a higher level of control over the time evolution of 
entanglement. In this case, even if the initial state is fixed, the 
variation of the coupling between the subsystems during time 
evolution may lead to many different paths for the entangle- 
ment to follow during the same interval of time. 

In this article, we propose the use of time-dependent Hamil- 
tonians to control the entanglement creation between qubits. 
We do it in an active way, i.e. after defining the shape wanted 
for the entanglement time evolution, the coupling magnitude 
is dynamically varied in order to achieve the expected behav- 
ior. The use of time-dependent Hamiltonians to establish en- 
tanglement is in fact gaining impetus in the last years. In 101, 
time-dependent spin-spin couplings are used for creation of 
long-distance entanglement where it is predicted an increase 
of entanglement between the first and last spin of a chain 
whenever the ac part of the coupling has a frequency match- 
ing the Zeeman splitting. In pO, sine and rectangular field 
frequency modulation is shown to be favorable to improve, 
enhance and stabilize the degree of the atom-field entangle- 
ment in a typical atom-cavity model. In 16], optimal control 
theory is used to determine the coupling modulation that leads 
to maximum logarithmic negativity for a pair of opposite os- 



cillators in a harmonic chain. We now provide a solution for 
the entanglement design between qubits through the choice 
of the appropriate wave form of the coupling intensity. This 
paper is organized as follows. In Section [III we provide an 
approximate solution for the inverse problem of obtaining the 
time dependent coupling from a desired time dependent en- 
tanglement function. In Section |III1 we discuss the meaning 
of these different paths and find applications when the qubits 
are subjected to decoherence. We summarize and conclude in 
SectionllVl 



II. STATEMENT AND SOLUTION OF THE PROBLEM 
A. Ansatz 

In order to illustrate the general idea behind entanglement 
design, let us consider the simple two-qubit XY Hamiltonian 

^W^^K^ + '^M), (1) 

where X{t) is a time-dependent coupling magnitude between 
resonant qubits and cr's are usual Pauli matrices. In the com- 
putational basis, if the initial state of the system is \ip{0)) = 
|01), the evolved state will be 

\iP{t)) = cos77(t)|01) -isin?7(i)|10), (2) 

where 

7?(t)= / A(r)dr (3) 

"'0 

is the area of the interaction pulse. Since the evolved state ^ 
is pure, the entanglement between the qubits can be quantified 
by means of the entropy of entanglement |[^. The result is the 
composite function 

sm] = -c[77(t)] \og,{c[rjm - log^isMm, (4) 

where c[77(<)] = cos^ ri{t) and s[77(i)] = sin^77(t). Our goal 
is to draw the time evolution of the entanglement function. 
Mathematically, we would like to obtain the area pulse rj{t) 
that satisfies the equation S[ri{t)] = f{t) for an arbitrary func- 
tion f{t). This function f{t) is the shape we want the entan- 
glement to have. For two qubits, this function is limited to 
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assume values between and 1, and for the initial preparation 
considered here, f{t) must also obey /(O) = 0. Should such 
r]{t) be found, an experimentalist able to control the interac- 
tion area pulse will not only gain control over the amount of 
entanglement created from t = to f = T, but also on how 
the entanglement actually evolved during this time interval. 

We provide now an approximate analytical solution for 
this problem, and illustrate its validity with some particular 
choices of f{t). Clearly, the pulse area will have to be a 
function of f{t), otherwise the problem does not make sense. 
Consequently, the entanglement itself might be considered a 
function of f{t). When thinking this way, the solution of 
S[f{t)] = f{t) is nothing more than an attempt to linearize 
the entropy of entanglement, regarding f{t) as the indepen- 
dent variable. This linearization (or approximate lineariza- 
tion) is precisely our strategy to solve the problem. 

We start by considering the linear entropy Sl ~ 2(1 — 
Trpp ) where pn is the reduced state obtained by tracing out 
one of qubits. Using we obtain 5*^ = sin^ [277(f)]. The 
design of the linear entropy can be performed exactly because 
the solution of Sl = f{t) is just 

77(t) = iarcsin[/(t)i/2]. (5) 

However, our main goal is to draw the entropy of entangle- 
ment. We then slightly modify the solution ^ trying to com- 
pensate the information lost when moving from the entropy of 
entanglement to the linear entropy. We propose the following 

ansatz 

77(t) = iarcsin[/(i)«/2], (6) 

where q is an adjustable parameter. We are now going to opti- 
mize the value of q such that S[i]{t)] becomes as close as pos- 
sible to f{t). We have already mentioned that this is equiv- 
alent to approximate yi{f;q) = S{f) to ?/2(/) = / when 
we think of / as an independent variable (/ <E [0, 1]). Let 
M = C[0, 1] the set of all continuous functions on [0, 1]. The 
functions j/i(/; q) and 1/2 (/) are elements of A/. We can de- 
fine a metric (distance) for this set which for j/i and 1/2 reads 

i 

d{q)^ f \yi{f;q)~y2{f)\df. (7) 
Jo 

In FiglT] this distance is plotted as a function of the adjustable 
parameter q. It is clear that there is a value of q{q ~ 1.345) 
that minimizes the distance between both curves to less than 
5 X lO^''. In order to analyze to which extent this is a good 
approximation, we show in Figl2]the entropy of entanglement 
as a function of f{t) using our ansatz (|6]l. An exact so- 
lution would lead to a perfect straight line with slope equal 
to 1 and crossing the origin (dashed line). One can see that 
our ansatz, with that value of q, is indeed a very good solu- 
tion of the problem. This plot shows that for small values of 
f{t), i.e. small entanglement, the solution works less well as 
indicated by the slight deviation from a perfect straight line. 
However, as we are going to see from explicit examples, this 
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Figure 1 : Distance between the entropy of entanglement and a per- 
fect straight line as a function of the adjustable parameter q. The 
inset shows the region where lies the minimum of d{q). 



deviation does not do much harm to the design of entangle- 
ment. The time dependence of the coupling is found from the 
first derivative of (|6|l with respect to time 
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Figure 2: The designed entropy of entanglement S as a composed 
function of the target function fit). The inset shows the region where 
the ansatz deviates more from a perfect straight line. 

In Figl3l we show an example of entanglement design 
where we choose the entanglement temporal shape to be 
f{i) = 1 — e""*. One can see that there is a very good agree- 
ment between S{i) (the designed entanglement) and fit) (the 
target entanglement). This is in accordance with the success- 
ful linearization shown in FiglS] However, this example brings 
out an important fact, namely the coupling constant diverges 
at times when the entanglement goes to zero (in this example, 
at t = 0). In fact, we can see that \{t) given by ^ will di- 
verge whenever f{t) = (because q/2 < 1) or /(<) = 1. 
These divergences seem to be consequences of approxima- 
tions. It is unlikely that the exact linearization (if possible to 
find) would lead to such divergences of the coupling for well- 
behaved functions. Besides, the form of the entropy of entan- 
glement and the unitarity of time evolution (causing ri{t) to be 
continuous) forbid the choice of discontinuous /(<). In gen- 
eral, physical instantaneous abrupt changes in entanglement 
result from measurements and not from unitary evolution. In 
spite of these divergences, the ansatz has been proved to offer 
a very accurate solution of the problem. In the next subsec- 
tion, we present a renormalized version of the ansatz which is 
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Figure 3: Example of temporal control of entanglement, (upper) 
The target function is plotted using dotted lines, and the designed 
entropy of entanglement is plotted using solid lines, (bottom) The 
time evolution of the coupling leading to the designed entanglement. 



Figure 4: Example of temporal control of entanglement obtained us- 
ing the renormalized ansatz. (upper) The target function is plotted 
using dotted lines, and the designed entropy of entanglement is plot- 
ted using solid lines, (bottom) The time evolution of the coupling 
leading to the designed entanglement. 



able to attain similar levels of precision but with the advantage 
of leading to physical (finite) couplings. 



B. Renormalization Algorithm 

As a first attempt to eliminate the divergences in A(t), one 
could think of considering q > 2. This would certainly re- 
move the divergence caused by /(O^^+^Z^ around f{t) = 0, 
but this would not help to solve the difficulties that arose at 
f{t) = 1. Besides, as seen from the distance function plot- 
ted in FiglT] the precision of the solution would be poor. A 
more effective idea is to keep q < 2, actually keep the opti- 
mum value q = 1.345, and then cut off the divergence of X{t) 
whenever it appears, i.e., around f{t) = and f{t) = 1. This 
lead to the simple renormalization algorithm 



m 
m 



1 qf{t)-'+^'^ 
Aq otherwise, 



df{t) 



dt 



if 5a<f{t)<5i 



(9) 



with Ao finite, 5q > 0, (5i < 1, and q = 1.345. FigH] shows 
the renormalized version of the previous example, with the 
choices Ao = 0, (5o = 10^^, and 5i = 1 — Sq. The parameters 
Ao, i5o, and 5i can be chosen as to minimize the error while 
keeping X{t) finite, but we leave optimization to be studied 
elsewhere. As one can see from FigH] there is still an excel- 
lent agreement between the target curve and the one obtained 
with our renormalized ansatz (|9]l. 

In Fig|5] we provide another example of active temporal 
control of the entanglement evolution. The target function is 
now the exotic f{t) = 5 + 7 arcsin[sin(7rKt — 7r/2)]. Again, 



there is a very good precision in the design of the entangle- 
ment. In both plots, we used q — 1.345, Ao — Q, 5q — 10^^, 
and 5i = 1 — 5o. 
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Figure 5: Example of temporal control of entanglement obtained us- 
ing the renormalized ansatz. (upper) The target function is plotted 
using dotted lines, and the designed entropy of entanglement is plot- 
ted using solid lines, (bottom) The time evolution of the coupling 
leading to the designed entanglement. 

Our intention in this paper is not to treat optimization prob- 
lems like minimization of the time taken for the system to 
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achieve a certain amount of entanglement. Instead, we show 
that the system may have its entanglement following a chosen 
path between two fixed values. The Hamiltonian considered 
in this work as well as initial system preparation fix the start- 
ing point to be ebit and allow the end point to be any value, 
including 1 ebit. 



C. Possible physical setting 

We would like now to analyze the feasibility of our idea 
in experimentally accessible physical systems. Experimental 
and theoretical work on superconducting qubits for quantum 
information purposes is a very fast growing field of investiga- 
tion ll9l- [Tlll . It is currently possible, for example, to experi- 
mentally couple two su perc onducting flux qubits through an 
Ising-type Hamiltonian il2l - [l4ll 

^ = - E 2^^^"^ + ^^''^^ + (10) 
1=1 

where and represent the bias and tunneling energies of 
qubit i, respectively, and J is the qubit-qubit coupling con- 
stant. 

There have been reported many different couplers between 
flux qubits providing a good control of the sign and magni- 
tude of the J coupHng M [TiHiell . In for example, it 
is theoretically shown how a monostable rf- or dc-SQUID can 
mediate coupling between two adjacent flux qubits described 
by Hamiltonian ( fTOl i. and the experimental implementation of 
this scheme is reported in ifTill . In this experiment, measure- 
ments demonstrate that J can be sign and magnitude tuned 
with a properly chosen flux bias in the monostable rf SQUID. 
This tunable coupling, added to the fact that flux waveforms 
can be electronically controlled, may be used to obtain an 
arbitrary time-dependent coupling J(t), the key ingredient 
needed in our scheme. In particular, the experiment reported 
in ifTsIl presents an improved tunable coupling element (com- 
pound Josephson junction rf-SQUID) able to provide a sign 
and magnitude tunable J{t) with minimal nonlinear crosstalk 
from the coupler tuning parameter into the qubits. 

For our purposes, Hamiltonian[TO]is as good as ([TJ as long 
as the tunelling energy is taken equal zero, and the qubits are 
initially prepared in the state |V'(0)) = |H — ), where |±) are 
eigenstates of the Pauli matrix ax- In this case, the interaction 
picture evolved state considering Hamiltonian ( fTOl ) will be 

\i^{t)) = cos7;(t)|+-) -ismTj{t)\-+), (11) 

which is locally equivalent to (|2]i. Consequently, the qubit- 
qubit entanglement is still given by ^ and our scheme can 
then be applied with no modifications. The case A 1.2 7^ 
may work too by including application of microwave drivings 



in. INTERPRETATION OF THE ENTANGLEMENT 
PATHS 

The question that naturally arises is: why would one 
want to obtain control over entanglement trajectories? This 
section is dedicated to a possible answer to this question. 
Let us consider the region defined by the ordained pairs 
{{Kt, Entanglement) £ [0,10] x [0, 1]}, just like consid- 
ered in the plots shown in Figs.© and (|5]l. The curve R{t) — 
{Kt/ 10) divides the plane into two different sets according to 
entanglement: trajectories above this curve mean that the sys- 
tem visits states with higher entanglement, while going from 
to 1 ebit, than when f{t) < R{t). In the ideal case, it is always 
possible to achieve 1 ebit, independently on the trajectory tak- 
ing place above or below the boundary defined by R{t). As 
we are going to show, different trajectories controlled by us- 
ing our ansatz provide a direct visualization of the effect of 
the decoherence acting on the qubits. This gives us a first clue 
about the meaning and application of different choices of X{t) 
(|9j for particular user-defined entanglement trajectories f{t) 
originally defined in the ideal case. 

In the previous section, we presented the qubit-qubit dy- 
namics following a designed coherent evolution to achieve 
a final state with maximum amount of entanglement. We 
now analyze the problem in the presence of two representative 
types of decoherence channels: the amplitude-damping (AD) 
and phase-damping (PD) channels. The description of the in- 
coherent dynamics of a composed system is usually made by 
using the appropriate master equations in the Lindblad form. 
The derivation of these equations are discussed with great de- 
tail in 1 17]. The density matrix p{t) describing the system 
state under Born-Markov approximations obeys the master 
equation in the Lindblad form 

m = -t[H{t),p{t)]+Y,[Lkpit)Ll-^LlLkp(t) 

k 

-]^p{t)LlLu], (12) 

where the first term gives the unitary part of the time evolution 
while the others describes the irreversible coherence loss into 
the environment. The operators Lk are the called Lindblad or 
quantum jump operators ifisll . 

The AD channel describes the dissipative dynamics where 
the qubit loses excitation to the environment, at decay rate F. 
Thus, each qubit dissipative process is described by a single 
Lindblad operator Lk = ^/2^(J^ , with cr^T = {af. — ia^)/2 
(fc = 1, 2). Substituting these jump operators in ( fT2] i. we have 
the system master equation under the AD channel 

p{t) = -i[H{t),p{t)]+ J2 n^cTkPitH -a+a^pit) 

fe=l,2 

-pitW^k]- (13) 

The PD channel is a purely quantum channel leading to 
losses of quantum coherence without any energy relaxation. 
In this channel, each qubit suffers the action of a single Lind- 
blad operator Lk = VFct^ (k = 1,2). By substituting these 
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jump operators in (fT2] i. we have the master equation for the 
PD channel 

m = -i[H{t),p{t)] + ^ V[alp{t)al-p{t)]. (14) 

fc=l,2 

To quantify the entanglement content of mixed states, we 
consider the entanglement of formation, which is defined as 

EoF ^h[{l + ^/l-C^)/2), (15) 

where h{x) = ~x\og2X — (1 — a;)log2(l — x) and C is 
the concurrence. For the initial preparation I V'(O)) — |01),the 
density operator obeying the master equations (fTSl l or ( fT4l i fall 
into the class of what is called X states lfl9ll . For these very 
symmetrical states, concurrence assumes the simple form 



C = 2max{0, |p23(t)| - TM^WW, ImWI 

-V/P22(i)p33(i)}, (16) 

where Pij{t), with = 1, ... ,4), are the corresponding 
matrix elements of p{i) in the computational basis. 

We are know in position to get back to the original question. 
In order to answer it, we will consider what happens with final 
amount of entanglement achieved under the influence of these 
decoherence channels, when different prescribed couplings 
\{t), originally obtained from user-defined functions f{t) 
ranging from to 1 ebit, are calculated with our ansatz 
We define a family of functions = (Kt/lO)^, where the 
parameter p selects the region delimited by R{t) = f{t; 1). If 
p > 1 [p < 1) the path in the ideal case would access states 
less (more) entangled before reaching 1 ebit, as delimited by 
R{t). In Fig. |6l we show the final amount of entanglement 
achieved when following a path given by f{t; p) under the in- 
fluence of each decoherence channel. From these plots one 
can see that the AD channel does not remove the degeneracy 
of symmetric paths p and 1/p. The situation changes radi- 
cally for the PD channel, where one can clearly see that it is 
more advantageous to spend most of the time in states with 
less entanglement. 

IV. FINAL REMARKS 



scheme for controlling \{t) may find applications in experi- 
ments aiming at entanglement maximization using XY Hamil- 
tonians in the presence of the PD channel. We hope this 
work can trigger studies including generalizations for more 
involved systems such as general spin chains or coupled con- 
tinuous variable systems. 



EoF 
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Figure 6: Final amount of entanglement achieved when following a 
path given by f{t) — in the presence of the AD channel 

(top) and PD channel (bottom), as a function of T/k (vertical axis) 
and (logj^Q p) (horizontal axis). 



To summarize, we have proposed a method to control the 
time evolution of entanglement between two qubits in real 
time. In order to achieve such control, we proposed the use 
of time-dependent coupling magnitudes, and provide a rela- 
tion between the interaction area pulse or coupling constant 
wave form and the desired shape of the entanglement func- 
tion. These entanglement trajectories, under specific condi- 
tions, revealed an independence of the final entanglement on 
each of the symmetric trajectories (defined in the text) for the 
case of amplitude damping. The phase damping channel has 
been shown to not possess this symmetry. Consequently, our 
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